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f  inal  Report  to  Rome  Am  Development  Renter  for  Sponsored  Research  m  "Oblique 
lonospnerte  Heating  ami  Effects  or,  Radio  Propagation"  during  Mav  i  !  S-ptember  .'JO. 
1987.  Prepared  by 


Dr.  Min-C’hang  Lee 
Principal  Investigator 


1.  Introduction 

Preliminary  research  under  the  sponsorship  of  the  RADC  has  been  actively  conducted 
at  Massachusetts  Institute  of  Technology  in  the  following  three  areas:  (1)  Self-focusing 
of  high  power  radio  waves  traversing  the  ionosphere.  (2)  Nonlinear  scattering  of  radio 
uo'o  off  the  ionospheric  density  striations.  and  (3)  Enhanced  ionospheric  heating  by  'he 
combined  operation  of  two  healers.  Briefly  described  as  follows  is  the  research  progress  in 
1  hese  projects. 

2.  Self- focusing  of  radio  waves 

A  graduate  student.  H.C.  Han.  has  analyzed  the  self-focusing  of  a  Gaussian  radio  wave 
beam  propagating  through  a  uniform,  magnetized  plasma  under  the  supervision  of  Dr. 
Min  ('hang  Lee.  This  work,  when  applied  to  the  ionospheric  radio  propagation,  offers 
the  following  scenario.  The  radio  waves  injected  from  a  ground-based  transmitter  can 
be  approximated  as  plane  waves  during  their  transionospheric  propagation.  Through  the 
filamentation  instability  (Lee  and  Kuo.  Radio  Science.  1985),  wave-like  ionospheric  density 
fluctuation-  can  be  excited  and.  consequently,  break  the  radio  wave  beam  into  filaments. 
These  filaments  ran  be  approximately  described  by  Gaussian  radio  wave  beams. 

The  -eif-focusing  process  under  consideration  arises  from  the  radio  wave- ionospheric 
plasma  interaction-  that  vioid  nonlinear  forces  acting  primarily  on  electrons.  These  forces 
mav  cause  the  increasing  of  ionospheric  plasma  dielectric  constant  along  the  radio  wave 
pat  li  and  results  m  the  focusing  of  radio  waves.  The  nonlinear  forces  induced  by  radio  waves 
are  Hie  thermal  pressure  force  and  the  ponderorriotive  force.  These  nonlinear  forces  work 


.iddit  ivelv  in  c  ausing  the  self- focusing  of  radio  waves  though  one’  may  become  dominant 


■  >\i'  Hie  other  under  certain  ionospheric  conditions.  For  instance,  the  thermal  pressure 
force  i-  t ire  dominant  force  when  the  scale  lengths  of  induced  plasma  density  fluctuations 
ah>ue  the  geomagnetic  field  line  exceeds  several  electron  mean  free  paths. 

under  the  seif-focusing  process,  the  shape  of  a  radio  wave  beam  can  be  altered.  More 
sp'Tiiicadv.  a  Gausdan  radio  wave  beam  will  change  its  field  intensity  profile  during  the 
•  ransionospheric  propagation.  The  enhanced  field  intensity  near  the  beam  center  due  to 
the  seif-focusir.g  process  can  result  in  intense  electron  healing  and  lead  to  the  generation 
nf  hot  spots  in  the  ionosphere.  As  the  first  step  in  our  theoretical  analysis  of  the  problem, 
we  oni\  consider  the  ponderomotive  force'  effect  for  simplicity  The  threshold  power  of 
the  concerned  process  is  determined  by  the  balance  between  the  natural  diffraction  and 
the  seif  focusing  process.  It  depends  upon  the  initial  beam  width,  the  nonlinear  dielectric 
constant,  and  the  wave  number. 

Solving  the  nonlinear  c  ave  equation,  we  can  clearly  me  t  fie  trend  of  the  spatial 
evolution  of  the  wave  field  intensity  profile,  namely,  a  Gaussian  radio  wave'  beam  doc- 
not  maintain  its  shape  during  trie  self- focusing  process  in  the  ionosphere.  If  we  ignore  the 
-  duration  pre-esses,  the  focal  length  can  be  defined  to  be  the  propagation  distance  before 
trie  wave  held  intensity  readies  infinity  due  to  the  self-focusing eifect.  This  theoretical  focal 
length  represents  the  upper  bound  of  the  characteristic  path  length  for  the  self- focusing 
process  to  accomplish.  Since  the  thermal  pressure  force  and  the  ponderomotive  force 
contribute  additiveiv  t< >  the  seT  foeusing  <>:  radio  waves,  our  calculated  threshold  power  is 
apparently  overestimated  because  iri  our  formulation  of  the  theory  we  ignore  the  effect  of 
the  thermal  pressure  force.  However,  interestingly  enough,  our  threshold  power  has  the 
same  order  of  magnitude  as  that  found  by  the  numerical  work  of  Bernhardt  and  Duncan 
I. VI  P.  HJ82j  who.  by  contrast,  neglect  the  effect  of  the  ponderomotive  force.  In  our  future 
iiunvtiiai  study  of  the  mif  focusing  process,  we  shall  include  the  thermal  pressure  force 
! o  obtain  a  complete  physical  picture  of  the  concerned  process.  Attached  with  the  report 
as  t  fie  Appendix  A  is  the  preprint  of  a  proceedings  paper  that  will  be  published  by  the 
Ionospheric  Effects  Symposium  held  during  5  7  May.  HJS7. 

3.  Nonlinear  scattering  of  radio  waves  in  the  ionosphere 

\  nonlinear  scattering  process  has  been  analyzed  by  Keit  h  M.  Groves,  another  graduate 
undent  of  Dr.  Min-('hnne  Lee  to  show  the  nonlinear  inode  conversion  of  radio  waves  into 
'•iei  i  rosi  at  ic  plasma  modes  in  the  presence  of  ionospheric  density  striations.  This  work 
v. as  motivated  by  the  observations  of  the  spectral  broadening  of  monochromatic  \  IT 
radio  ognals  traversing  the  ionosphere  conducted  in  the  so-called  YLF  wave  injection 
experiments  (Bell  et  ai..  1083). 
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Wcoruing  lo  Bel!  el  id.  \  1983).  a  bandwidth  expansion  as  large  as  100  Hz  was 
experienced  bv  a  nearly  monochromatic  signal  at  13. (i  KHz  :  i  Hz  injected  from  a  ground- 
based  YLF  transmitter  as  the  signal  traversed  the  ionosphere  and  reached  the  satellite 
;-b  nudes  iri  the  range  of  bOl)  3S0O  kilometers.  In  other  words,  the  expansion  oi  bandwidth 
a?  large  as  I'd  of  the  carrier  frequency  was  observed  in  the  transionospheric  propagation  of 
YLF  waves.  I  he  otT-carrier  components  arc  believed  to  be  electrostatic  in  nature  (Inan  and 
Bell,  i 9S5 ) .  This  phenomenon  occurs  only  in  the  presence  of  impulsive  YLF  hiss  and/or  a 
lower  hvbrid  resonance  (LHR)  noise  band  with  an  irregular  cutoff  frequency,  and  only  for 
dgnals  whose  frequencies  exceed  the  LHR  frequence-  ;it  the  satellite  location. 

Two  mechanisms  have  been  proposed  to  explain  the  spectral  broadening  phenomenon. 

One  is  the  nonlinear  scattering  process  of  YLF  signals  by  the  existing  ionospheric  density 
fluctuations  that  renders  the  nonlinear  mode  conversion  of  YLF  waves  into  lower  hybrid 
v.avi  -  The  scattering  of  YLF'  waves  by  ionospheric  density  fluctuations  causes  elliptica.ll v 
:>oi arizen  modes.  However,  we  find  that  the  induced  ellipticaily  polarized  modes  may 
become  predominantly  electrostatic.  They  occur  when  the  scale  lengths  of  ionospheric 
density  irregularities  are  much  less  than  the  YLF  wavelength.  \\>  note  that  the  scattered 
YLF  waves,  therefore,  have  a  broadened  wave  vector  spectrum  that  can  give  rise  to  the 
Doppier  shift  of  the  YLF  waves  as  the  movin'-  satellite  sensed. 

The  other  proposed  mechanism  suggests  parametric  instabilities  that  ran  produce  a 
spectrally  broadened  signal.  The  injected  YLF  wave  excites  both  the  Stokes  and  anti- 
Stokes  components  of  lower  hybrid  waves,  concomitantly,  producing  a  field-aligned  purely 
growing  mode  Lee  and  Kuo.  ;9S4  .  The  Doppler  shift  of  the  sidebands  relative  to  the 
whistler  pump  wave  created  a  spectrally  broadened  signal.  A  paper  submitted  bv  the 
graduate  student.  Keith  M.  Groves  and  accepted  already  for  the  Student  Prize  Paper 
( '•  .;n petition  in  the  1988  National  Radio  Science  Meeting.  Boulder.  Colorado  is  attached 
as  the  Appendix  B  of  this  report. 

4.  Enhanced  ionospheric  heating  by  two  heaters 

The  other  project  conducted  by  I)r.  Min-Chang  Lee  at  MIT  with  the  participation 
of  graduate  students.  C.P.  Liao.  D.R.  Rivas,  and  K.M.  droves  is  the  theoretical  analysis 

,,f  enhanced  ionospheric  heating  bv  the  combined  operation  of  two  heaters.  I  he  proposed  - - 

>r 

,  i ■ ; i a i  io  i--  a.s  Mllows.  Illuminate  the  ionospheric  F  region  bv  a  vertically  injected  BF  or  - 

B 

ML  boater  wave  first.  Then,  launch  a  YLF  wave  propagating  through  the  HF  or  MF  rt 

modified  ionospheric  region. 

it.  ims  been  known  both  theoretically  and  experimentally  that  short-scale  (typically 
"■ricr  Mnlrj  ionospheric  density  striations  can  be  excited  bv  a  vertically  injected  11F  or 

/ 
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MI'  wave  from  a  ground-based  transmitter  in  less  than  a  second  (Lee  et  a!..  1986).  while 
larce-'i  ait  i  on  .  I  in  ml  rods  of  meter  scale)  striatums  need  tens  of  seconds  to  he  excited. 

1*  is  desirable  to  operate  the  UK  or  MF  heater  in  a  pulsed  mode  with  a  duration  of. 

-a\.  .!  lev*,  seconds  to  assure  the  favorable  excitation  oi  meter-scale  rather  than  large- 
scale  ionospheric  density  striations.  \s  discussed  in  the  previous  section,  the  HI'  or  MF 
wave  induced  meter-scale  ionospheric  density  striations  are  able  to  noniinearly  scatter  the  » 

subsequently  launched  YLF  wave  into  a  quasi-electrostatic  mode.  viz.,  a  lower  hybrid  wave 
mode.  « 

Several  prominent  ionospheric  effects  can  be  expected  and  tested  experimentally.  The 
most  notable  effects  cue  associated  with  the  generation  of  lower  hybrid  wave  modes  that 
can  effectively  accelerate  electrons  along  the  geomagnetic  field.  The  electron  acceleration 
lead-  to  airglow  (6300  A.  5577  .1.-  •  ■)  enhancement  and.  then,  a  broad  height  distribution 
of  plasma  lines  that  can  be  measured  by  incoherent  backscatter  radars.  Presented  in 
A  upend’ x  C  i-  trie  research  progress  in  this  work. 

1  he  further  theoretical  investigation  of  this  project  shall  be  pursued  and  extended 
to  t’ne  case  that  the  HF  waves  are  injected  obliquely  from  a  ground-based  transmitter. 

Collaboration  with  L)r.  Stanford  Yukon  of  RAU(‘  is  expected.  Dr.  Yukon's  ionospheric 
radio  research  group  has  been  actively  planning  an  oblique  ionospheric  heating  experiment 
with  other  Air  Force  research  groups  and  universities.  During  the  oblique  ionospheric 
heating  bv  an  HF  wave,  the  short-scale  ionospheric  density  striations  cannot  be  excited 
because  the  radio  wave-ionospheric  plasma  resonance  conditions  cannot  be  matched  for 
the  excitation  of  parametric  decay  instability.  However,  relatively  large  scale  (>  tens 
of  meters)  ionospheric  density  striations  can  he  generated  via  other  plasma  instabilities 
-uch  a-  filamentation  instability/ self- focusing  instability.  Brillouin  instability,  and  Raman 
instability.  Hence,  in  principle,  these  instabilities  are  still  capable  of  converting  the  VLF 
waves  into  lower  hybrid  waves  via  the  nonlinear  scattering  process. 

The  calculation  of  wave  field  intensities  in  the  caustic  region  of  an  obliquely  propagating 
HF  .cave  were  carried  out  by  Field  and  Warber  (1985).  1  heir  calculations  show  that  a 
transmitter  having  a  power-gain  product  of  5  MW  can  launch  an  oblique  wave  strong 
enough  to  produce  electric  fields  of  several  tenths  of  a  volt  per  meter  nr  more  in  the 
caustic  region.  If  the  wave  field  intensities  are  indeed  so  intense  as  Field  and  Warber  * 

i  19X5  i  indicate,  strong  ohmic  heating  of  electrons  can  be  expected.  Nevertheless,  in  the 
recent  oblique  experiment  conducted  in  Maine  by  AFC1L.  RADC.  1  Diversity  of  Lowell, 
ami  .MIT.  no  indication  of  electron  heating  was  seen.  \\i  shall  look  into  this  problem  in 
the  near  future  for  the  oblique  heating  experiments  currently  being  planned  by  the  RADC 
and  A  F  ( II.. 
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"i  I  resrntation  of  research  results  in  conferences 

I  'An  i  rim  made  fur  at t ending  o  >nierenees  were  supported  oy  the  R  A  1 )( ’  c<  >n  tract .  One 
.'.a-  the  trip  t<>  Springfield.  Virginia  in  May.  POST  made  by  Dr.  Min  Chang  Lee  with  four 
-indents  in  present  five  papers  in  the  Ionospheric  Kffects  Symposium.  I  lie  other  trip  was 
made  bv  1  Jr.  MinChanu  Lee  alone  to  \  ancon ver.  Canada  for  delivering  an  invited  tall;  at 
•  he  1 A  Cl  A  (ieneral  Assemble  :n  August.  !  ‘AST . 
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YP  PEND  IX  A 


THE  SELF  FOCUSING  OF  A  GAUSSIAN  RADIO  WAVE  BEAM 


A  bstract 


\s  a  high- power  radio  wave  beam  propagates 
!i  the  I'.niwphere.  the  nonlinear  interactions  between 
'tie  wave  and  plasmas  can  lead  to  the  self-focusine  ut 
•he  radio  beam.  The  required  lioeshold  or  critical 
power  is  determined  by  ’he  balance  between  the  natural 
diffraction  and  the  self-focusing  process.  Solving  the 
nonlinear  wave  equation  near  the  wave  beam  center,  we 
can  obtain  the  focal  length  for  the  concerned  process  to 
accomplish.  It  has  been  shown  that  a  Gaussian  radio 
wave  beam  shall  not  maintain  its  shape  during  the  self- 
focusing  process  in  the  ionosphere. 


1.  Introduction 


I  lie  seif  focusing  of  radio  waves  m  'pare  has  re¬ 
ceived  much  attention  m  the  research  on  ionospheric 
modification  bv  high  power  radio  waves.  This  phe¬ 
nomenon.  has  also  been  actively  studied  in  optical  me¬ 
dia  see.  e  g.  Chiao  et  al..  1964;  Kelley.  1965;  Lu- 
govoi  and  Prokhorov.  1974  and  references  therein  .  In 
ionospheric  plasmas,  the  self- focusing  of  rad’o  waves 
arises  from  the  wave-ionosphere  interaction  that  yields 
nonlinear  forces  acting  primarily  on  electrons.  These 
forces  may  cause  the  increasing  of  ionospheric  plasma 
dielectric  constant  along  the  radio  wave  path  and  re- 
-ults  in  the  focusing  of  radio  waves.  The  nonlinear 
forces  induced  by  radio  waves  are  the  ponderomoti ve 
force  and  the  thermal  pressure  force.  When  the  scale 
length  of  plasma  density  perturbation  along  the  geo¬ 
magnetic  field  line  exceeds  several  electron  mean  free 


App-A-1 


ail.-. 


the  'i.ermai  pressure  1 1 >r <  f *  dominates  nwr  the 
;>  ail  ;<!'■ 'in.  a  i  \  <•  tone  n,  causing  the  'fit  focusing  nisii 
1  ' !  i  1  t  \ 

l  n>if  r  t  li<-  'f i i - locusi nu  p:oce>'.  i  hr  shape  of  a 
.niii.  ''f.tm  fin  be  altered.  More  'pedt  naily . 

,t  t.aii'siaii  radio  beam  wiii  change  it'  field  intensitt 
proiii*'  tiirouu.ii  t  tie  t  ransionosphenc  propagation,  [lie 
purpose  "i  the  present  work  is  to  investigate  analvtu  aliv 
'It'  prohlem.  A'  the  first  step.  we  restrict  ourselves 
on  tjje  [ion  ileromot  l  ve  force  eftect  onlv.  One  can  see 
•  he  trend  ot  the  spatial  evolution  of  the  field  intensity 
profile.  The  threshold  fieid  for  the  self-focusing  process 
’tirtis  out  to  he  dependent  of  the  initial  beam  width, 
the  nonlinear  dieted  ru  i  oust  ant.  and  the  wave  riuinher. 
From  the  analysis  of  the  intensitv  protile.  certain 
‘."iinipt tons  niaiie  m  solvinu  the  wave  equation,  such  a- 
'he  expansion  ot  the  held  intensity  as  a  function  of 
near  the  beam  miter,  and  the  neulect  of  the  diffraction 
ft  ted  .  ■  an  ■  fust  it  ied 


2.  Nonlinear  Wave-Ionosphere  Interactions 

When  a  i.'juh  power  radio  wave  Ffam  propagate' 
a.  t  :.»•  !"ti"- phere .  ’he  radiation  pressure  utves  rise  to 
a  tr  aiiitif  ar  1>  «rce.  Known  as  the  pomif  romot  i  ve  force, 
.ft  i ; .  u  hasicall’.  on  electrons,  if  the  wave  frcquenct 
mm  ti  htulier  than  the  ion  plasma  frequency.  1  he 
;>•  >ndf  ri-mot  ive  torce  acting  •  m  an  electron  is  given  b\ 

f  v  /  '  ^  1 :  <  1  - 

Im.d 

where  .1  i'  tlie  amplitude  of  *iie  wave  electric  field. 
/•.  .  1(  r  ■  i  o'  jji.  'nice  thus  force  points  in  t  lie  tlireci ton 

■  ip p< is 1 1 e  to  that  of  VI'.  it  i'  able  u>  push  the  elect r< ais 
out  of  1  tie  wave  pat  ll . 

\i:  fb-drou  under  the  action  of  ponderomot i ve 
! 1 1 r *  f  t f i ■  i -  an  apparent  fleet  rtcai  potential  determined 


I  mid  ’he  i  •  unlit  ions  that  1'  :s  dowiv  varying  in  time 
and  :t  ha-  a  -•  ah-  ff.iuth  much  ureater  than.  t||f  Debvi 


'■n'j’i..  .  •  i  that  the  >  !..<  ail  \  •  i.e-mali/ed.  the 

Hi  -It  ,m:. ,iiin  relation  in  tin'  presence  nonr.eromot  1  ve 
n  >rce  '  "m  "■  '.vrit  ten  a-* 
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a  iit'ii'  •  .  tilt-  liHcksirotnui  electron  density.  k  p  :>  tin- 

H<  >it  znwuin  constant.  F,  is  nm  temperature,  and  /,  i- 
eie.-t  n<;*  • ' uiperat  ure.  Iberejore.  the  dielectric  .r.slant 
'  a  ’he  plasma  is  niven  by 
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a  lie  r<’  t  ;  -  tii'‘  i  iielec  t  nc  <  •  >U>  t  a  lit  ter  ! : Tf  -  Dace .  ,t  tl  (i 
tiie  at' r.  uri  mini  elec  t  ti  iti  plasma  lretiuen<’\.  \  -sumniu 
'iii'  Aa'.r  treiiuency  _  !-  iimr:i  l, -eater  than  •  :,<■  incai 
e|ei p.a-ma  iretpiencv  we  can  express  t'.q.  I1 
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15.  .N oitiinear  Wave  Equation 


'  •  •  1 1  - 1 1 i e r  a  powertm  radio  wave  propauatnm  in  / 
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[f  we  assume  that  i ( j  *  )  •  ■ we  orih  need 

to  r «■  i ;ti 1 1  tit*-  irst  nvo  terms  m  Eq.  ■  ■'!.  1  lien  ncglec*  ng 
t  he  iirst  term  in  Eq.  ilOt  v;e|d> 
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l.ti  /  t'\  r.  t  i  exp  nkz  w  1  )  .  where  i  i~  the  >  1  •  •  v\  1  \  - 

varying  field  amplitude  vec.or.  k  is  t  he 

unperturbed  wave  number  ami  u.1  is  tin  \va\e  frequency, 
riien.  we  i  ave 
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Il.e  quant  it ie>  A  .  A  .  T  and  -  are  m-fined. 
respectively.  as  the  sca.c  length  oi  t  aiong  ti.e  i  axis, 
that  perpendicular  to  zaxis.  t  lie  wave  cugtli.  the  wave 
period,  and  'he  time  scale  of  the  temporal  variation  oi  t  . 
f  or  simplicity,  we  consider  a  cylindrical  wave  beam  with 
T  ■■  ~  which  implies  the  quasi  stationary  condition, 

and  \  A  --  A  .  The  order  of  magnitude  estimates 
of  each  term  in  Eq.  (12)  are  made  as  follows; 
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Neglecting  the  higher  order  terms  but  keeping  the 
transversal  part  of  the  first,  the  second  and  the  third 
tenm  in  Eq.  I  12 )  leads  to 
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4.  Profile  Evolution  of  Wave  Field  Intensity 


l.et  i'  \  1.  o.  :  i  exp  / _-|  o.  :  i  Substituting  t  into 
Hu.  i  1 i  ami  separating  the  real  ami  imaginary  parts  oi 
i  l.t..  we  have  the  following  two  equations 
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expressed  m  t  lie  cylindrical  coordinate  system  with 
azimuthal  symmetry  for  a  uniform,  isotropic  collisionai 
plasma.  Ihese  two  equations  can  he  rewrit.en  as 
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Eq.  ilTi  can  be  further  rewritten  in  the  form  of 
the  Hamilton- Jacobi  equation  e>r  'he  particle  motion: 

-21’  =  0.  where  ;  ( T5'r  -  21  is 
the  Hamiltonian.  V S  is  the  momentum  and  1  is  the 
potential  field  Wagner  et  al.  1 9 1> .  We  can  specify  the 
time  with  kz.  the  spatial  coordinates  with  ko.  >  with 

.  and  \  with  -  .1  s—  —  .  The  first  term  with  \ 

“  ‘  k *  A  '  p 

represents  the  diffraction  effect  and  the  second  term 
represents  the  nonlinear  effect  li.e.  the  focusing  effect  i. 
from  the  Hamilton-.] acobi  equation,  we  have 
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This  etjuation  describes  the  radial  position  o  of  a  ray  as 
a  junction  <4  ■  >n  it-  trajectory,  as  an  analogue  ’<>  the 

particle  motion  m  a  potential  field.  The  trend  for  the 
spatial  evolution  of  the  wave  field  intensity  profile  can 
‘  hen  t  >e  seen 


Suppose  that  we  have  a  (lausstan  heam  with  t  he 
f:,eid  inteu'itv  tin)  Tmrx:>i  1  .  u  where  lTn  -,s 


the  peax  intensity,  and  u\,  is  the  beam  width.  Since 
f  '  i_>.4~(expt  1 I  Irom  hq.  it'i.  the  potent  tai  l  is. 

lienee,  a  junction  ot  n  given  bv 
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1  he  iirst  two  terms  on  the  right  hand  side  represent 
diffraction  while  the  third  term  represents  the  focn.-mg 
effect,  if  l’  has  a  maximum  at  o  -  gmax-  then  t he 
focusing  effect  is  dominant  for  g  <  gmaj.  whereas  the 
diffraction  effect  is  dominant  for  p  ■  g„iar-  1  he  p,„<ij 
can  be  determined  by  taking  the  derivative  of  \  with 
respect  to  g  as  follows 
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'p 


(21 ) 


[lie  "critical  power  for  self- focusing  can  be  defined  by 
'Citing  g,nnJ  —  0  in  Eq.  (21  I.  viz.. 


(221 


Consider  a  Gaussian  beam  with  .4m  7$>  -4m  and 
a  plane  wave  front  (i.e.  ^  ;=o  =  0).  The  portion 

of  the  beam  at  p  <  gmax  will  be  focused,  while  that 
at  p  •  Omar  will  be  diffracted.  In  other  words,  an 
immediate  redistribution  of  field  intensity  is  expected. 
The  shape  of  the  beam  will  change  from  a  Gaussian  to 
one  with  a  sharper  center  and  higher  wings. 

In  fact,  for  any  field  profile  peakine  at  the  renter 
and  dying  out  as  p  gets  large,  we  can  do  the  local 
analysis  near  n  --  0*.  We  can  expand  the  field  intensity 
profile  in  a  series  of  p  for  p  —  0~.  Then,  this  scries  ran 
be  written  as  1  -  qn  ~  •  •  •  after  normalizing  the  first 
term  and  defining  a  new  variable  p.  Substituting  this 
series  into  Eq.  (19),  we  can  see  that  n  can  not  be  less 
than  2  to  avoid  the  physically  meaningless  singularity  of 
T  at  g  —  0.  With  this  restriction,  together  with  a  high 
enough  beam  power,  the  focusing  effect  will  dominate 
over  diffraction  near  the  beam  center,  Ender  the  afore¬ 
said  conditions,  it  is  then  possible  to  solve  Eq.  (13) 
analytically  near  the  beam  center. 
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5.  Self- Focusing  Near  The  Beam  Center 


lor  simple  mat  liematical  manipulation,  we  express 
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where  .S'lp.ci  =  5r<  p.  z  i  -  tS^g.zi.  If  the  wave  power 
is  much  greater  than  the  critical  power,  we  can  make 
the  following  approximations: 

i  i  t  The  diffraction  terms  in  Eq.  ( 2d )  can  be  neglected 
as  o  —  I) 

( ii  l  Sr{  o.  : )  —  unf  ;  i  -  g~aM  -  )  -  r  I  -  •  •  • 

iiii)  By  assuming  a  spherical  wave  front,  we  have 
\ i  p.  z  >  ■  bf)l  ~  \  -  g'b-jl  z ) 

I’sing  the  above  expansions  for  Sr  and  S,  in 
Eqs.  i2dl.  (24)  and  collecting  terms  of  the  same  order 
in  p.  we  then  get  a  set  of  equations  for  o' s  and  6s. 
If  we  neglect  terms  higher  than  the  second  order  of  n. 
Eqs.  r2di.  and  (24!  become 
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<J 6,  -  ^-6oao  —  0 

( 28  i 

where  arn  -  e~a°.  If  we  define  a  function  g{ z)  = 
J~  bn  z')<lz'.  then  Eqs.  |26)-(2S)  have  the  forms  of 
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where  <  anti  <  _>  art-  integration  constants.  I  he  initial 

COIltil  t  loll'  <t  If  set  to  lie 
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wtiere  R  is  the  initial  curvature  of  wave  front  at  ;  =  0. 
and  R  ■  0  for  a  diverging  beam,  while  R  -  0  for  a 

converging  beam,  flere.  u\,  is  the  initial  beam  width 
and  l, n,j,  is  the  initial  peak  field  intensity.  Then.  Eq. 
i  26*  I  can  be  straightforwardly  solved.  In  the  following, 
we  express  our  final  results  in  terms  of  the  squared  peak 
field  internin'  and  the  scjuared  beam  width  as  functions 
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Setting  / 1  a  i  --  0.  we  can  determine  the  focal  length  as 


Amo  h  2 

ll'n  \  ( j,  R 
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where  tp  and  .ij  have  been  defined  in  Eqs.  (61  and 
(9l.  rcspectivel v.  These  results  were  obtained  by 
Akhmanov  et  al.  1966.  for  studying  the  self-focusing 
effect  in  nonlinear  optical  medium.  But  the  assumption 
of  a  global  quadratic  field  intensity  profile  and  the 
geometrical  optical  approximation  are  required  in  their 
work.  From  >>ur  above  analysis,  one  can  see  that  both 
of  those  are  unrealistic  and  unnecessary  for  the  near¬ 
axis  solution.  For  example,  in  the  ionospheric  F-region. 
n  l.o  ■  10iJm  T,  '■  F,  --  1000°K.  and  the  wave 
frequency  is  taken  to  be  100  MHz.  wn  ~  100m.  then  the 
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critical  iield  intensity  is  found  tn  he  .F.('  ■  3d."Y  m.  For 

a  beam  with  .4m  —  10.4',/,.  the  local  length  calculated 
from  Fq.  i  34 )  is  about  4.2  Km. 


6.  Summary  and  Discussion 


In  summary,  our  analysis  has  showed  that  a 
Gaussian  beam  no  longer  maintains  its  shape  during 
the  self-focusing  process.  A  threshold  or  critical  power 
is  required  for  the  self- focusing  effect  to  overcome  the 
natural  diffraction.  If  the  incident  power  is  much 
greater  than  the  critical  power,  the  beam  will  focus  at 
the  beam  center.  The  theoretically  determined  focal 
length  represents  the  upper  bound  of  the  characteristic 
path  length  for  the  self- focusing  process  to  accomplish. 
This  is  because  we  have  not  considered  any  energy 
absorption  process  that  may  impede  or  cause  the 
saturation  of  the  concerned  self-focusing  process.  The 
fact  that  the  phase  and  amplitude  of  the  wave  field 
vary  drastically  near  the  focal  point  as  shown  in 
Eqs.  >25)  and  (32)  also  shed  doubt  on  the  validity  of 
neglecting  the  higher-order  terms  in  Eq.  (12).  The 
thermal  effect,  expected  to  contribute  additivelv  with 
the  ponderomotive  force  to  the  self-focusing  effect, 
should  also  be  taken  into  account.  More  detailed 
analysis  of  this  problem  will  be  presented  later. 
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APPENDIX  3 


SPECTRAL  BROADENING  OF  VLF  RADIO  SIGNALS 
TRAVERSING  THE  IONOSPHERE 


K.M.  Groves 

Research  Laboratory  of  Electronics 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 


Abstract 


Nearly  monochromatic  signals  at  13.6  kHz  del  Hz  injected  from  a  ground-based  VLF 
transmitter  can  experience  a  bandwidth  expansion  as  high  as  1%  (  100  Hz)  of  the  incident 

wave  frequency  as  they  traverse  the  ionosphere  and  reach  satellite  altitudes  in  the  ranne 
of  600-3800  kilometers  Bell  et  al. ,1983].  We  investigate  two  different  source  mechanisms 
that  can  potentially  result  in  the  observed  spectral  broadening  of  injected  monochromatic 
VLF  waves.  One  is  the  nonlinear  scattering  of  VLF  signals  by  induced  ionospheric  density 
fluctuations  that  renders  the  nonlinear  mode  conversion  of  VLF  waves  into  lower  hybrid 
waves  These  quasi-electrostatic  modes  result  when  the  injected  VLF  waves  are  scattered 
by  ionospheric  density  fluctuations  with  scale  lengths  less  than  0. 7( c/  fp){  fe  /  f0 ) 1 ' ' .  where 
c .  fp.  ft,  and  Jo  are  the  speed  of  light  in  vacuum,  the  plasma  frequency,  the  electron 
cyclotron  frequency,  and  the  VLF  wave  frequency,  respectively.  A  second  mechanism 
involves  the  excitation  of  electrostatic  waves  (low'er  hybrid  waves,  low  frequency  quasi¬ 
modes  I  by  the  injected  VLF  waves.  This  process  tends  to  produce  a  spectrally  broadened 
transmitted  pulse  with  peaks  at  a  discrete  set  of  frequencies  on  both  sides  of  the  nominal 
carrier  -quenev  Controlled  study  of  the  spectral  broadening  phenomenon  bv  the 
combined  operation  of  two  lone  HF  and  one  VLF)  heaters  in  the  so-called  ionospheric 
modification  experiments  is  suggested. 

Thesis  Advisor:  Dr.  M.  C.  Lee 
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1.  INTRODUCTION 


In  this  paper  we  propose  two  mechanisms  to  explain  the  recently  observed  bandwidth 
expansion  experienced  by  nearly  monochromatic  signals  at  13.6  kHz  :r  1  Hz  injected  from  a 
ground-based  VLF  transmitter  as  they  traverse  the  ionosphere  and  reach  satellite  altitudes 
in  the  range  of  600-3300  kilometers  (Figure  i)  Bell  el  al..  1 983 j -  A  schematic  illustration 
of  this  puzzling  phenomenon  is  shown  in  Figure  2. 


FIGURE  1.  Geometry  of  spectral  broadening  observations. 


FIG  I  RE  2.  Bandwidth  expansion  of  transionospheric  VLF  wave 


This  expansion  of  bandwidth,  which  resuits  in  a  proportional  increase  in  signal-to- 
noise  ratio,  may  be  as  large  as  i'T  ( —  100  Hz)  of  the  carrier  frequency,  and  the  on-carrier 
components  are  thought  to  be  eiectrostatic  in  nature  Inan  and  Bell,  1 9S5  .  The  five 
tvpes  of  spectrally  broadened  transmitter  signals  which  have  been  found  are  illustrated  in 
Figure  3.  They  occur  only  in  the  presence  of  impulsive  YLF  hiss  and,  or  a  lower  hybrid 
resonance  i  LHR)  noise  band  with  an  irregular  cutoff  frequency,  and  only  for  signals  whose 
frequencies  exceed  the  LHR  frequency  at  the  satellite  location.  The  observed  signal  at  the 
carrier  frequency  is  usually  suppressed  relative  to  the  sidebands. 

The  iinear  scattering  source  mechanism  was  first  suggested  by  Bell  et  a/.,  .1983!,  who 
hypothesized  the  creation  of  the  required  ionospheric  density  fluctuations  by  precipitating, 
low  enerev  {<  1  kev  i  electrons.  Such  precipitation  events  have  been  accompanied  by  both 
YLF  mss  and  irregular  LHR  noise  bands  (McEwen  and  Barrington,  1967;  Laaspere  et  al.. 
1971;  Gurnett  and  Frank,  19721.  Bell  et  al.  (1983)  then  speculate  that  the  bioadening 
of  the  transmitted  pulse  spectrum  results  from  the  scattering  of  the  initial  signals  from 
the  precipitation-induced  density  fluctuations  and  the  subsequent  coupling  into  quasi- 
electrostatic  whistler  mode  waves  of  short  wavelength.  The  Doppler  shift  associated  with 
these  short  wavelength  modes  is  large  enough  to  produce  the  bandwidth  expansion  of  the 
signals  measured  on  a  moving  satellite. 

The  first  mechanism  which  we  propose,  presented  in  section  2,  differs  from  that  of  Bell 
et  al.'s  in  two  basic  aspects:  1)  The  density  irregularities  required  for  scattering  may  be 
induced  by  the  injected  VLF  waves  themselves,  rather  than  particle  precipitation  :  Lee  and 
A'uo.1984);  2)  The  proposed  scattering  is  the  nonlinear  scattering  of  VLF  signals  by  induced 
ionospheric  density  fluctuations  that  renders  the  nonlinear  mode  conversion  of  VLF  waves 
mto  lower  hybrid  waves.  In  other  words,  the  scattering  of  VLF  waves  by  ionospheric 
density  fluctuations  causes  eiliptically  polarized  modes.  The  induced  elliptic a.11  y  polarized 
modes  mav  be  predominantly  eiectrostatic. 

The  second  mechanism,  outlined  in  section  3,  produces  a  spectrally  broadened  signal 


signal  level  (40)  signal  level  (40) 


via  parametric  instabilities  in  two  steps.  In  the  first  step,  the  injected  VLF  wave  excites 
both  the  Stokes  and  anti-Stokes  components  of  lower  hybrid  waves,  concomrrutantlv 
producing  a  field-aligned  pureiv  growing  mode.  The  lower  hybrid  waves,  acting  as  pump 
waves  .  can  then  interact  with  low  frequency  thermal-level  density  fluctuations,  generating 
both  bacKScattere--  ■  nd  forwardscattered  lower  hybrid  wave  sidebands.  The  frequency  shift 
of  the  sidebands  tive  to  the  lower  hybrid  pump  wave  creates  a  spectrally  broadened 
signal.  The  processes  may  be  summarized  as  follows: 

Step  1  i  VLF  - •  LH]-z  —  Zero  Frequency  FA  Mode 

Step  2  i  LH)  - -  LHn~  -  Low  Frequency  FA  Mode 

where  VLF.  LH.  ana  FA  stand  for  Wry  Low  Frequency  wave.  Lower  Hybrid  wave,  and  Field 
Aligned,  respectively.  The  -  -  indicates  that  both  upshifted  and  downshifted  modes  are 
produced  The  excitation  of  LH  waves  by  these  processes  may  aiso  result  in  the  acceleration 
of  b ~'f h  electrons  and  ions.  A  summary  and  conclusion  follow  in  section  4. 

2.  NONLINEAR  WAVE  SCATTERING 

Theory 

A  monochromatic  \  LF  wave  transmitted  from  a  ground-based  station  into  space  has 
been  observed  to  change  from  linear  into  circular  polarization  ii.e..  whistler  mode)  as  it 
Taveis  through  the  neutral  atmosphere  into  the  ionosphere  Kintner  et  al  .  ’.983  If  a 
ducted  whistler  wave  mode  is  considered  for  simplicity,  the  wave  electric  field  may  be 
represented  as 

E0  =  E0(  x  w  iy  i  exp  if  k0z  -  wnt  I  1  I  ■ 

where  ’he  :  axjs  has  been  ’  aken  along  the  geomagnetic  field  lines:  is  the  transmitted 

\  ave  frequency  and  kn  the  associated  wave  vector,  assumed  to  be  alone  the  :  axis  for  a 
* ur ’ e-i  mode,  as  ilbist rated  in  Figure  la.  Propagating  into  an  unperturbed  ionosphere 


'on -0-  j 


the  whistler  wave  satisfies  the  electromagnetic  wave  equation 


'•'■here  jG 
lonospher: 
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allowing  : 

where  ajs 
W  ith  ■  h  e 
.  rreeuiar: 
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is  the  uniform  oscillatory  current  driven  by  the  incident  wave  held  in  the 
c  plasma. 

i  e "  N  o  E  0 

Jo  =  lb  1 

m,tw0  -  <4) 

m,  and  fi.  are  the  electron  charge,  uniform  background  plasma  density,  electron 
the  unsigned  electron  cyclotron  frequency,  respectively. 
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•!(;i  RF.  1  (leometrv  of  wave  propagation  for  ai  uniform  ionosphere; 
hi  ionosphere  with  irregularities  present 


presence  of  field -aligned  ionospheric  density  irregularities,  the  scattered  wave. 
>r  longitudinal  modes,  can  be  described  bv 
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-the  .uniform  current  induced  by  the  interaction  of  the  whistler  wave  field 
densitv  irregularities  We  take  the  x  axis  along  the  direction  of  the  density 
tics  which  are  assumed  to  be  of  the  form 


An  -  Afi  exp  i k x 


depicted  graphically  in  Figure  4b,  The  nonuniform  current.  6j»,  can  then  be  expressed  as 

<\)s  =  -et  N0cv  -  — j0j  i  5i 

•’  O 

where  fv  is  the  induced  velocity  perturbation.  The  scattered  wave  field  has  the  general 

form  of 

Es  =  xE»  -  iy  E_v .  exp  n  kc2  —  h  i  6 1 

Solving  the  wave  equation  together  with  the  electron  momentum  equation,  the  scattered 
wave  field  is  found  to  be  eilipticallv  polarized 


Es  -  E0(  (''n  N0 )  x  -  ( x  -  iy  )ES|  exp  it  kcz  -  a,v>t  h 


F-s  is  the  amplitude  ot  the  circuiarlv  polarized  component  of  the  wave  given  by 

•->  ,2  . 
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where  vp,  and  c  are  the  electron  plasma  frequency  and  the  speed  of  light  in  vacuum, 
respectively,  and  *  ;s  the  wavenumber  of  the  density  perturbation  defined  by  (41.  Writing 
the  scattered  field  in  terms  of  the  ionospheric  irregularity  scale  length.  A  =  2rr  k,  We  find 


where. 


—  1 

A c  ~  v  2t(c,  i( o )  ^  4.4  k0 


The  scattered  held  is  composed  of  a  linearly  polarized  i  LPl  part  and  a  circuiarlv  polarized 
<  P  part  as  indicated  in  >9).  If  the  ratio  (A*  A")  «<  1  ii.e..  k  ':>>  k0  i  the  <catter<"2 
wave  ;s  dominated  by  the  linearly  polarized  component,  which  oscillates  m  the  direction 
•  l  ’he  irregularity  wave  vector,  k;  hence,  the  wave  is  electrostatic  in  nature  and  exhibits 
i  nroanenecl.  "mareed  wave  vector  spectrum  relative  to  the  incident  wave.  The  popp.er 
•hift  frequency  measured  bv  a  moving  sn'e'I'te  will  be  enlaced  as  vveil. 


I  •inniiii 


The  condition  found  above  for  effective  nonlinear  scattering  into  quasi-eiectrostatic 
modes  requires  \  <<  Ar.  Typical  plasma  parameters  for  the  upper  ionosphere  lower 
magnetosphere  region  of  interest  are  ^^,2: r  =  0.3  MHz.  FD  2tr  =  0.6  MHz.  Assumine 
an  incident  frequency  nt  v,,  2~  —  13.6  kHz.  we  calculate  Ac  a-  1.3  km.  The  magnitude  of 
'Me  bandwidth  expansion  measured  by  a  moving  satellite  due  to  the  Doppler  shift  of  the 
-•'attered  VLF  wave  is 
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v’nere  vs  is  the  velocity  of  a  satellite  moving  across  the  earth’s  magnetic  field  as  shown 


in 


Figure  5.  and  ks  -  kx  -  k0z. 


FIGURE  5.  Schematic  of  Doppler  shift  phenomenon 

Given  v,  ^  »  km  sec  and  A/  ^  100  Hz.  we  find  A  rs  SOm  <<  \c.  so  the 
..ndition  for  nonlinear  scattering  is  well  satisfied.  Thus  we  have  found  that  ionospheric 
.-regularities  with  scale  lengths  of  several  tens  of  meters  ran  produce  Doppler  shifts  of 


•he  observed  magnitude.  Ionospheric  irregularities  can  occur  naturally  or  be  created  by 
•he  incident  VLF  waves.  Particle  precipitation  can  produce  irregularities  with  a  broad 
range  of  scale  lengths  in  a  region  near  the  F-laver  peak,  about  300  km  above  the  earth's 
-urface.  However,  the  reduced  phase  speed  of  the  scattered  quasi-elertostatic  wave  leads 
•o  sigiuncant  Landau  damping  of  this  wave  mode,  and  the  amplitude  of  such  a  wave 
decreases  exponentially  as  it  propagates  upwards  to  satellite  altitudes  of  600-3800  km. 
Nevertheless,  the  required  irregularities  can  be  induced  by  intense  VLF  waves  locally  Lee 
and  Kuo.  1984  .  Furthermore,  t  hese  wave-induced  irregularities  possess  a  relatively  narrow 
spectrum  of  scale  lengths  ranging  from  a  few  meters  to  a  few  tens  of  meters. 

3.  EXCITATION  OF  ELECTROSTATIC  MODES 


Parametric  Excitation  of  Lower  Hybrid  U’ares 


In  this  section  we  discuss  the  proposed  second  mechanism  responsible  for  the  observed 
VLF  spectral  broadening  that  involves  the  parametric  excitation  of  lower  hybrid  waves  and 
ionospheric  irregularities.  As  in  the  previously  considered  case  of  nonlinear  scattering,  the 
incident  wave  is  also  assumed  to  be  a  ducted  whistier  mode  with  k0  =  k0z.  propagating 
along  the  geomagnetic  field  lines,  taken  to  be  the  :  axis.  Assuming  the  space-time 
dependence  of  the  perturbations  to  be  exp  i(kc  •  r  -  wtL  ,  the  first  physical  process  of 

VLF  - LHj-  -  Zero  Frequency  FA  Mode  can  be  described  by  the  following  wave 

frequency  and  wave  vector  matching  relations: 


O  —  \ 


kj.  -  ks  =  zk0  =  kj_  -  k,  (111 

.-.  *  he  --ubsrr-.nts  -  and  s  refer  to  the  Stokes  anti-Stokes  components  of  the  high 

freouencv  iower  hybrid  sidebands  and  the  field-  aligned  zero  frequency  mode,  respectively, 
'’’hooding  ’he  r  axis  to  coincide  with  the  wave  vector.  k«.  of  the  field- aligned  modes  .  as 


-down  :n  Figure  Ha.  and  expressing  the  frequency  w,  —  m  {Re'w, 


=  01.  the  matching 


conditions  i  11  i  may  be  represented  by 


-v  i  x  —  -"o  ttt 

ki  _  =  zk0  —  xks  i  1 2) 

This  process  has  been  investigated  in  detail  in  Lee  and  Kuo.  1984;.  For  the  upper 
ionosphere  in  the  frequency  range  of  interest,  the  nonoscillatory  current  resulting  from 
the  beating  of  the  density  fluctuations  associated  with  the  lower  hybrid  waves  and  the 
electron  response  to  the  whistler  pump  field  is  found  to  be  the  dominant  nonlinear  effect. 
The  scale  lengths  of  the  zero  frequency  modes  are  found  to  be  a  few  tens  of  meters;  hence, 
its  kQ. 


FIGURE  6.  Illustration  of  wave  vector  matching  relations  for  a)  Step  1;  b  1  Step  2 


Stimulated  Scattering  of  Lower  Hybrid  fUares 

The  lower  hvnd  waves  generated  by  the  first  process,  with  wave  vector  and  wave 
frequency  denoted  bv  kj.  -j\.  respectively,  can  act  as  pump  waves  to  produce  high 
r e ( i ! i e n c v  lower  hvbrid  sidebands  and  low  frequency  field-aligned  modes  ik,.-,. 


.\od-H-  1  0 


Re  W,  <<  _c,  i  via  the  second  parametric  instability  described  by  the  following  matching 
relations. 

-*■'  2  ±  =  1  —  -*-’j 

\<2~  =  kj  r  ks  1  13) 

where  the  subscripts  -  -  refer,  to  the  anti-Stokes ,  Stokes  components  of  the  lower  hybrid 
sidebands.  The  wave  vector  matching  relation  for  this  step  is  illustrated  graphically  in 
Figure  6b.  The  Stokes  component  of  the  sidebands  may  be  chosen  as  an  eigenmode  of 
'he  plasma  by  setting  ks  =  2k*,  which  minimizes  the  threshold  of  the  instability.  While 
the  backscattered  sideband  mode  has  the  same  scale  length  as  the  pump  wave,  that  of 
'he  forward  scattered  sideband  mode  is  one-third  of  the  pump  wavelength.  The  dispersion 
relation  obeyed  by  the  LH  pump  wave  is  given  by, 
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is  the  lower  hybrid  resonance  frequency. 

The  high-frequency  sideband  modes  are  driven  by  the  beating  current  induced  by  the 
pump  wave  field  on  the  density  perturbation  associated  wilh  the  low-frequency  mode.  The 
nonosciilatory  low- frequency  mode,  by  contrast,  is  driven  primarily  by  ion  nonlinearity, 
including  the  ponderomotive  force  and  the  nonosciilatory  beating  current.  While  the 
ponderomotive  force  acts  as  the  driving  source,  the  nonosciilatory  beating  current  imposes 
a  stabilization  effect  on  the  instability. 

The  plasma  dynamics  under  consideration  are  determined  by  the  momentum  equation. 


mjUji  d '  d\  -  vj  ■  VWj  =  mq^E  -  vj  *  B0)  -  VTjnj 


•  14) 


■fie  >";  i<r.ion  oi  continuity. 


i  <9/  <’dt  )m  -  V  -  njVj  =  0 
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and.  in  the  present  case  of  low  frequency  oscillations,  the  quasi-neutrality  condition, 

nj  =  n>  i 16) 

where  n.  q,  and  T  represent  particle  density,  charge,  and  temperature  in  units  of  energy. 
The  subscript  j  refers  to  the  particle  species:  i  for  H~  ions  or  e  for  electrons,  as  written 
explicitly  in  (16).  Wavelike  perturbations  are  assumed,  having  the  form  <fn].V]  -- 

expukj  r  -  aiM  :  n , ,  v,  exp  t(ks  ■  r  -  ^j,t )’;  and  8n~ ,  v*  exp  t(k2~  •  r  -  -;2*  t ) 
where  k2±  and  are  defined  by  (13).  Under  the  electrostatic  approximation  we  may 
also  express  the  electric  field,  E  .  as  the  gradient  of  a  potential  function.  <z>.  The  interaction 
of  the  hieh  frequency  LH  pump  wave  with  the  high  frequency  LH  sidebands  drives  the  low 
frequency  field-aligned  modes:  the  corresponding  momentum  equation  is. 

vsmjvsj  =  iqj/n0(vsj  <  B0 )  -  ksTj(nS/ no )  -  i Fj  (17) 

x  here  Fj.  the  nonlinear  ponderomotive  force  arising  from  the  convective,  has  the  form 

Fj  -  -  irnij  k0  ivVvoj  -  v  _jV  *  j )  -  L-vVv,j  -  k_  vojvV 

Appivmg  (16)  to  the  equation  of  continuity,  we  obtain 

w5(ns;  n0 )  =  k,  •  (vse  -  v"1)  =  k5  i  vsi  -  v;nl)  (  IS) 

where  the  induced  electron  and  ion  nonoscillatory  beating  currents  are  defined  by 

novj11  =  8nmOJv+  j  -  k/b  -  Sn^jV^  -  Sn'_}voj  (19) 

The  velocity  responses  of  the  electrons  and  ions  to  the  lower  hybrid  wave  fields  in  (  19)  may 
be  found  bv  solving  the  linearized  momentum  equations  for  both  species;  the  result  is 

veo  ^  =  — iko,±  x  X  -i-  iy(ui0,±/ne)  —  iz(ko*rie/ k0 ^-Uo.t ), 

, 

<  (  eOo,*  ■  mPnf  !  -  (  v“e<fneo,±  ,  nonf)j  -O' 

v;0  ^  =  'k0.±  *  x ( Q , /  ^o. -  —  ty  —  ii(k0i/  k0i  )) 

X  (e<P0.-  Tn,u;0,t )  -  (v“6nio  ^/no^o.j. )!  '  -1 
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The  equation  of  continuity  can  then  be  linearized  and  solved  for  the  corresponding  density 
perturbations. 
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'  (1  -k Lxv2/^±)-l  ,23) 

where  <50  =  0  and  6*  =  I.  Substituting  these  first-order  results  and  solving  (17W19) 
together  with  a  similar  set  of  equations  governing  the  interaction  of  the  LH  pump  wave 
and  the  low  frequency  field-aligned  mode,  and  assuming  w,  «C  u;Q,  the  following  dispersion 
relation  is  obtained 

a^5  —  but,  -  r^s  —  d  =  0  <  — -4  ) 


where 
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Substituting  u.,  =  -  i*y  into  (24).  two  equations  for  the  growth  rate.  .  and  the  real 

frequency,  ujr,  are  obtained. 

■v "  =  d,  (b  -  2aujr)l  -  -jj~  (26) 

i  b  —  2axjr  i3  —  b(  b  —  2au;r )“  -  ac(  b  ■+•  2au.’r )  -  a"d  =  0  (  27  i 

Setting  ■>  =  0  to  determine  the  threshold  of  the  instability,  (26)  and  (27)  may  be  solved 
for  the  coefficient,  d.  giving 

d  =  2a{  i  b  3a r  -  c ,  3a,  (  b ,  3a —  c  ' 3 a i 1 '  2  -  (b/3a)3  ~  be  6a‘  }  (28) 


Equating  this  result  with  the  expression  in  (25),  we  find  the  threshold  condition  of  the 
instability  requires 

vii  S:  v?hl  -  {}(ksxC5/^1)(2ki/ks)(  1  ~ 

<  (1  *  h)3'2  -  ^ h 1  “( 1  -  2h/3)]/ ;3(  1  -  ks2Cs2,fl2);5}C2  <  ;7?/a  -  -29. 

w  here  v,,  -  eE] ,  mja/j  is  the  ion  quiver  velocity,  h  =  b"  ;3ac,  and  Cs  =  (Tf  -  T,)  m,:  *  is 

the  acoustic  speed  in  the  ionospheric  plasma.  When  the  lower  hybrid  pump  field  strength 
is  sufficient  to  produce  quiver  velocities  in  excess  of  that  required  by  (29),  positive  growth 
rates  for  the  instability  occur. 


Characteristics  of  the  Instability 

This  mechanism  has  been  analyzed  for  typical  plasma  parameters  in  the  upper 
ionosphere  lower  maenetosphere.  namely  for:  -^pt :  2t  =  0.65A/  H  z\  fu  2t  -  0.55  MHz. 
Vp,  ,.2t  -=  1 5  KHz-.  Q,/2tr  =  :>30  H  =  \  =  5.45  *  10“4:  T*  =  T,  =  1000°K.  The 

analysis  was  carried  out  in  ’he  spectral  range  about  ks  2k0,  consistent  with  previous 
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assumptions.  Illustrated  in  Figure  T  is  the  spectral  dependence  of  both  the  growth  rate 
of  the  instability  and  reai  frequency  of  the  iow  frequency,  field-aligned  mode.  The  real 
frequency  of  this  iow  frequency  mode  corresponds  to  the  detected  broadening  of  the  VLF 
wave  frequency  by  the  orbiting  satellite  as  shown  by  the  frequency  matching  condition. 

13).  The  plots  in  Figure  7  were  calculated  by  solving  the  complete  dispersion  relation.  (221 
and  I  23).  for  the  instability  assuming  an  electron  quiver  velocity.  v0e  ;r  voi(^0/  f>, )  =  C\ 
The  growth  rate.  goes  to  zero  near  k,  =  2k0(  1  —  <).  because  the  threshold  electron 
quiver  velocity  exceeds  Cs  at  this  point.  As  can  be  seen  in  Figure  6.  ’he  growth  time  of 
the  instability  is  about  100  ms:  for  pulse  lengths  of  the  order  c;  one  second,  significant 
growth  can  occur.  The  growth  rate  and  spectral  range  of  the  instability  increase  linearly 
with  the  pump  power 

as  vJt  >  C's  v  2.  Although  the  calculated  magnitude  of  the  broadening  is  less  than  the 
observed  values  of  100  Hz.  it  increases  by  a  factor  of  two  for  an  order  of  magnitude  increase 
in  the  electric  field  strength  of  the  lower  hybrid  pump  wave. 

4.  SUMMARY  AND  CONCLUSION 

In  -ummary.  we  have  investigated  two  possible  source  mechanisms  that  can  be 
responsible  for  'he  observed  spectral  broadening  of  injected  VLF  waves.  In  the  presence 
of  mnospheric  irregularities  with  scale  lengths  of  several  tens  of  meters,  the  nonlinear 
■  cattering  of  the  VLF  waves  off  these  density  irregularities  can  produce  quasi-electrost atic 
modes  with  larger  wave  vectors  which  give  rise  to  the  appaient  spectral  broadening  through 
’he  Doppler  shift  observed  by  the  moving  satellite:  the  broadening  produced  through  this 
mechanism  la  of  'he  -ame  magnitude  as  the  observed  values.  The  amplitude  of  the  quasi- 
-iertrost  atm  modes  depends  linearly  on  the  amplitude  of  the  ionospheric  irregularities.  One 
possible  source  of  the  field-aligned  ionospheric  irregularities  is  particle  precipitation  Bell. 
,  !  ui..  . qs?,  ,  However,  irregularities  with  the  required  scale  lengths  can  also  be  induced 

-it-  1  n 
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locally  by  the  VLF  waves  themselves  via  the  parametric  instability  described  by  Lee  and 
Kuo.  1984  . 

Spectral  broadening  of  the  incident  wave  packet  may  also  stem  from  the  second 
mechanism,  which  involves  the  parametric  excitation  of  lower  hybrid  waves  and  low 
frequency  quasi-modes.  The  frequency  shift  resulting  from  this  -ocess  is  generally  less 
than  the  observations.  However,  it  is  significant  and  it  increases  with  the  field  intensity  of 
t  he  lower  hybrid  pump  wave.  These  two  suggested  mechanisms  can  contribute  additivelv  to 
the  observed  spectral  broadening  of  injected  VLF  waves  reported  by  Bell,  et  a!.,  ;  1 983  and 
others.  The  broadening  introduced  by  this  mechanism  is  significant  in  that  it  produces  a 
wave  frequency  spectrum  which  exhibits  suppressed  field  intensity  at  tne  carrier  frequency 
and  enhanced  intensity  at  a  discrete  set  of  frequencies  on  either  side  of  the  carrier  frequency; 
such  spectra  are  frequently  observed.  Furthermore,  the  lower  hybrid  waves  excited  by  this 
mechanism  have  enlarged  wave  vectors,  so  apparent  broadening  due  to  the  Doppler  effect 
comparable  to  that  calculated  in  the  first  mechanism  is  expected. 

Another  interesting  phenomenon  associated  with  the  excited  LH  waves  is  the 
acceleration  of  both  ions  and  electrons  in  the  ionosphere.  The  LH  pump  wave  generated 
in  the  first  step  of  the  mechanism  has  a  relatively  large  phase  velocity  and  can  interact 
effectively  only  with  electrons  in  the  tail  of  the  velocity  distribution  function.  By  contrast, 
the  phase  velocity  of  the  forwardscattered  LH  sideband  wave  generated  in  the  second 
step  of  the  mech  anism  is  one-third  of  the  pump  wave  phase  velocity,  and  the  parallel 
i  perpendicular )  component  is  within  an  order  of  magnitude  of  the  electron  (ion)  thermal 
velocity.  Thus,  the  upshifted  LH  wave  may  accelerate  bulk  electrons  and  ions  along  and 
across  the  geomagnetic  field  lines,  respectively. 

The  proposed  mechanisms  can  be  tested  in  the  so-called  ionospheric  modification 
experiments  with  the  <.  ombined  operation  of  two  heaters  in  the  following  scenario,  depicted 
:i  Figure  Initially,  a  high  frequency  t  HF)  heater  wave  is  transmitted  from  the  ground 
’n  illuminate  the  ionosphere,  producing  chort  ( meter  l-scale  field-aligned  ionospheric 
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FIGURE  8.  Ionospheric  Modification  Experiment 


irregularities.  The  subsequent  injection  of  a  VLF  wave  through  the  HF  wave-mod>fied 
ionospheric  region  wiil  induce  lower  hybrid  wave  modes  via  either  noniinear  scattering 
off  ionospheric  irregularities  or  parametric  instabilities,  as  elucidated  in  the  previous 
sections.  The  induced  lower  hybrid  waves  can  effectively  accelerate  electrons  and  thereby 
heat  the  ionospheric  plasma.  Physical  observables  associated  with  this  ueating  in'-'nd* 
airglow  enhancement  (e.g.,  6300  A,  5577  A.  and  shorter  wavelengths)  and  a  broad  height 
distribution  of  enhanced  plasma  lines.  In  addition,  backscatter  radar  measurements  can 
be  used  to  monitor  the  expected  intensification  of  the  field-aligned  ionospheric  density 
irregularities  by  :  ne  \  LF  waves. 
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IONOSPHERIC  MODIFICATIONS  BY  TWO  HEATER  WAVES 


Abstract 

The  ionospheric  modification  caused  by  an  HF  or  MF  heater  wave  can  be  enhanced  with  the  subsequent 
illumination  of  the  ionosphere  by  a  powerful  VLF  wave.  The  proposed  scenario  of  ionospheric  modifications  by 
the  two  heater  waves  is  based  upon  the  following  physical  processes.  Let  the  HF  or  MF  heateT  be  operated  in 
a  pulse- wave  mode  to  assure  the  excitation  of  short-  rather  than  large-scale  ionospheric  density  irregularities. 
These  excited  ionospheric  density  stnations  can  effectively  scatter  the  VLF  wave  into  a  lower  hybrid  wave 
via  the  nonlinear  mode  conversion  provided  that  the  scale  lengths  of  lonospnenc  irregularities  are  much  less 
'han  the  wavelength  of  the  VLF  wave.  For  example,  the  wavelength  of  a  VLF  wave  at  the  frequency  of 
.0  kHz  is  of  the  order  of  500  meters  in  the  ionospheric  F  region.  The  preferential  excitation  of  meter-scale 
'.nosDnenc  irregularities  bv  the  HF  or  N1F  heater  wave  can  provide  the  subsequently  injected  V  LF  wave  with  a 
favoraoie  condition  for  the  nonimear  mode  conversion.  These  density  stnations.  in  fact,  can  also  be  intensified 
by  the  powerful  VLF  wave  via  a  plasma  instability  that  can  concomitantly  generate  lower  hybnd  waves. 
The  ionosphere  modified  by  the  two  heater  waves  is  expected  to  have  intense  lower  hybnd  waves  and  short- 
icale  lonospnenc  density  stnations.  These  VLF  wave-produced  electrostatic  waves  can  effectively  heat  the 
lonospnenc  plasma.  Enhanced  modification  effects  in,  for  instance,  airgiow  and  height  distnbution  of  plasma 
lines  are  expected.  The  proposed  expenment  can  provide  the  controlled  study  of  the  spectral  broadening  effect 
y  propagating  V  LF  waves. 


1.  Introduction 

discussed  m  (.1  roves  e>  a!  ’.'»*•?  at  the  l *.>*» 7  Ionospheric  Effects  Symposium.  VLF  waves  ran  be 
r.. it. earn  scattered  .r.to  w*r r  nno  waves  ov  snort -scale  otiospnenc  density  irreguianties.  V  LF  waves,  it 
,  twnse  enough.  can  aiso  excite  lower  hybrid  w  aves  and  ionospheric  irregularities  via  parametric  instabilities 
’  e»  and  Kuo.  i‘lJS-1  Based  upon  'nesr  theoretical  wonts,  we  propose  an  experimental  scheme  tor  ionospheric 
'.rating  two  neater  waves  One  is  an  HF  or  MF  neater  wave  and  the  other  one  is  a  V  LF  heater  wave. 
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The  oropos^d  scneme  is  sj  loilows:  illumination  of  the  lonosDhere  by  an  HF  or  MF  wave  first  and. 
*hen,  a  VLF  wave  suDaeouentiy  i  see  Figure  i I:  has  oeen  mown  Doth  theoretically  and  experimentally  mat 
snort-scaie  u.e..  less  than  a  few  meters  i  lonospnenc  irregularities  can  be  excited  within  a  second  by  HF  heater 
waves,  wtuie  it  takes  tens  of  seconds  or  ionger  for  large-sc&ie  i  say,  hundreds  of  meters  and  longer  i  lonosDCenc 
.rreguianf.es  to  be  generated  see.  e.g.,  Fejer.  1979  and  Gurevich.  1978  for  a  review  .  As  shown  in  Groves  et 
al.  1987  and  elaborated  in  Section  2  of  tins  paper,  short-scaie  lonospnenc  density  stnations  can  etfectiveiv 
ause  the  nonlinear  scattering  of  VLF  waves  into  lower  hybrid  waves.  Hence,  we  oniy  need  to  operate  the 
HF  or  N1F  heater  m  pvusea  moae  to  assure  the  excitation  of  meter-scale  ionospheric  regularities.  During 
he  vertical  lonospnenc  neatmg,  the  HF  heater  wave  frequencies  should  be  less  than  the  F0F2  for  overaense 
heating  of  the  lonospnenc  r  region.  When  an  MF  heater  is  used,  the  wave  frequency  is  required  to  maun 
the  local  electron  gyrolrequency  in  tne  ionosphere  to  produce  short-scaie  lonospbenc  irregularities  via  electron 
-yciotron  seating  Lee  et  ai..  1986.. 

The  proposed  expenment  has  duai  purposes:  1  I  i  enhanced  lonospbenc  modifications  by  the  combined 
; Deration  of  two  neater  waves,  ana  I  2 )  controlled  study  of  the  spectral  broadening  ol  VLF  waves.  The  suggested 
diagnoses  tor  the  expected  lonospnenc  effects  shall  be  discussed  in  Section  3.  In  pnncipai.  the  two  mecnamsms 
-lucidated  :n  Groves  rt  al.  1987  can  be  ex  pen  mentally  distinguished. 


2.  Heater  Waves-Indueed  fonospheric  Effects 


\  anous  oiasma  instabilities  ran  be  excited  by  HF  heater  waves  dunr.g  overdense  ionospheric  -/Mina.  In 
mar.  a  ^ec.-nd.  =  r.ort -scaie  ‘Oospneric  density  striattons  can  be  generated  concomitant. v  won  Langmuir 
.Dper  ::vbn<3  waves  etc  requires  tens  of  seconds  for  the  generation  of  large-scaje  lonospnenc 

rregMiant-.es  by  seif- focusing  instability  or  tnet^i*'  fil»mentation  instability  Large-scale  irregularities  •  an 
give  rise  to  pnase  and  amplitude  scintillation  of  beacon  satellite  signals.  These  wave  interference  patterns  rr.av 
r.ot  be  associated  with  significant  attenuation  of  the  radio  signals.  By  contrast,  me  short-scaie  irr^guian  e*  are 
ante  to  cause  anomalous  absorption  of  the  radio  signals  via  nonlinear  scattering  that  converts  electromagnetic 
wav*  energy  into  electrostatic  wave  energy  (eventually,  plasma  kinetic  energy  ,  in  the  ionosphere 

3asicaiiv.  snort  scae  inRosonertc  irregularities  induced  by  an  HF  or  MF  heater  wave  can  scatter  me 
ircuiariv  poianzea  \  !.F  wave  into  an  elliptic  ail  y  polarized  wave  A  '.meariy  poianzed  component  d  tne 
altered  wave  is  .ntroduced  ~y  me  heid-aiigned  ionospheric  irregularities.  In  other  words.  me  scattered 
wnistie:  \  LF  wave  rieid  ;*  '-'-imposed  of  two  parts: 

Et  =  [‘Z  -  '-U'E-p  -  y£,p} 

E  ?  iiid  E  r  --nresent  me  wave  held  intensities  ■>!  the  •*ircu.,ari>  poianzed  component  and  me  .sr.ear.v 
>,*:■.  i  +  <\  om:>--d.*nt .  respectively  For  simplicity,  a  ducted  whistler  wave  propagating  a*ong  tne  ;  tion 

.as  ■  ***  b  -•  and  E  v  are  ‘oumi  to  be  r*:a!ed  ov 


*.  .ir*  i.*  •  •:  .**t>r»erir  r  r*gu.«ir. » :e«;  m-t  •  r.e  wa-,*;*r.gt  n  •  •  t  *  r.e  *  *  \\  + 

••«:  .  I:  a  ■  I.r  w  *-.e  •  <.  me  !r*«:*ier.-  '■  -l  i  i  *  Hz  is  m  !e«',*a.  and  *  lie  ■  *  nos  on  eric  •-••mlif  *  -d 

- e  *.  e.eriron  ''vr:otron  frequency  -  .}  MHz  ana  G  electron  p.asma  frequency  n  MHz.  the  wnisiier 

wave .engm  m  •  ::  ■  *»  d  n  e  r  *  -  '-..und  to  be  about  LOO  rr.e'ers  According  to  Equation  2  .  £~,p  .iomicates  .  vr 
.  ;  -  -  ;  *  -  at  ■  !!*:-•  -  meter-*.  aie  irr*g':.ar;t;es  produced  bv  HF  >f  MF  heater  wave  \:\  v 

--n<;e:  •  v  •;  .<:*  '  •.  *rs:on  *»f  a  v  LF  wnist  ier  '  wave  into  a  .j  wer  nvbrni  wave  I  he  :>ara-.e:  ar.d 

:  -rfe-.d.  ..  ir  *  a .  e  "  '  -s  #f  •  wer  n  v  o  r:  d  wave  are  identical  to  tne  wave  vector  ■  i  the  ducted  wmsi  .e* 
a\v e  •. . :  *  ;,.it  :  m*  ^e.,i  A  . g-e  *  .r.< -sDnertc  .rr^g  narities.  respectivriv 


.  -  ear  •>.  it  *  e ; g  ^roce<<  ..'■'■•.irs  at  ar.v  held  intensity  ol  the  incident  V  l.F  wave  -.Mramr*. r*-‘ 
e  ggere.-j  -  r.ter.se  enough  \l.F  waves  as  suggested  in  Lee  .ana  Kuo  and  (if  v  *s 

•■r*  -  -  air  •  nosnr.eric  ..’regularities  ivs  >nated  with  zero- frequency  modes  can  :-e  et.-  * e«; 

we-  - a  waves  -  v  •  “e  V  {  J-  \  aves  »  he  excited  i-vwer  hy  brid  waves  g:  w 

-i  -  ustat/iif.v  that  generates  daugnter  ower  hybrid  waves  and  low- ireouen-  v 


•  e  a 


:  ar  ari.e* 


ionospheric  irregularitie*  of  the  order  of  a  few  hert*.  Therefore,  the  HF  or  MF  heater  wave-induced  ihort- 
scale  ionospheric  irregularities  can  be  further  strengthened  by  the  V'LF  heater  ware*.  It  should  be  mentioned 
that  tero- frequency  modes  behave  as  a  standing  wave  pattern,  whereas  low-frequency  field-aligned  modes  can 
propagate  across  the  earth's  magnetic  field. 


Figure  1.  The  proposed  experiment:  illumination  of  the 
ionosphere  by  an  HF  or  MF  wave  first  and.  then, 
a  V'LF  wave  subsequently 


3.  Diagnosis  of  Expected  Ionospheric  Effects 

It  is  generally  believed  that  spread  F  echos  in  lonograms  and  scintillation  phenomena  are  caused  by  large- 
scale  (say,  Hundreds  of  meters)  ionospheric  irregularities.  The  generation  of  large-scale  ionospheric  density 
irregularities,  therefore,  can  be  sensed  by  ionosondes  and  scintillation  measurements  of  radio  stars  or  beacon 
satellite  signals.  By  contrast,  the  presence  of  field-aligned  short-scale  (a  few  meters  and  lessl  irregularities  can 
be  detected  by  backsratter  radars. 

3acitscauer  radar  measurements  can  also  be  made  for  detecting  the  excited  lower  hybrid  waves  when 
the  radars  beam  angle  is  j0°  with  respect  to  the  earth's  mignetic  field.  The  recorded  incoherent  backscatter 
'adar  spectrum  should  look  like  a  double-humped  ion  spectrum  peaking  at  the  lower  hybrid  wave  frequency. 
Since  lower  hybrid  waves  can  accelerate  electrons  effectively  along  the  geomagnetic  field,  airglow  enhancem  ^t 
due  to  the  imoact  excitation  of  neutrals  bv  energetic  electrons  can  be  expected.  Measurements  of  ajrgjow 
at  n.300  .4  l  reoi.  ~~  .4  '  green  i.  and  even  shorter  wavelengths  may  possibly  be  observed.  In  addition,  the 
'ecnnique  described  'n  Carlson  et  al.  1 982  for  finding  the  height  distribution  of  enhanced  plasma  lines  is 
capable  of  distinguishing  further  the  electron  acceleration  by  Langmuir  waves  and  lower  hybrid  waves  that  are 
produced,  respectively,  by  HF  and  VLF  heater  waves. 
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•4.  Comments  on  the  Proposed  Experiments 


There  are  several  ionospheric  seating  faculties  operated  at  the  frequencies  of  a  lew  MHz  in  tne  HF  or  N1F 
band  in  tne  I  S. A..  Europe,  and  the  I'.S.S.R.  VLF  transmitters  have  been  used  at  different  locations  tor  the 
study  of  wave-particle  interactions  in  the  magnetospnere  However,  it  may  not  be  possible  to  use  the  existing 
HF  'or  MF  i  and  VLF  transmitters  for  the  proposed  experiments.  Since  VLF  waves  in  a  narrow  frequency 
oand  can  be  generated  by  lightning  storms,  the  experiments  can  be  earned  out  with  available  lonosppenc 
heating  facilities  during  lightning  storms.  It  is  not  necessary  to  have  lightning  storms  occurnng  nearby.  If 
lightning  occurs  at  the  conjugate  location  in  the  opposite  hemisphere,  the  lightning-induced  VLF  waves  can 
propagate  in  ducted  whistler  wave  modes  and  bounce  back  and  forth  for  several  cycles  before  they  die  out. 

For  the  controlled  study  of  the  spectral  broadening  of  V’LF  waves,  the  expenments  require  an  air-borne 
\  LF  wave  detector  earned  by  a  moving  satellite  or  a  rocket,  flying  through  the  ionosphere  that  has  been 
modified  by  an  HF  or  MF  heater  and  a  VLF  heater.  The  predicted  spectral  broadening  effect  wiil  be  most 
prominent  if  the  aircraft  moves  across  the  geomagnetic  field  perpendicularly. 
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RADC  plans  and  executes  research .  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
C ommumcations  and  Intelligence  (C3T)  activities.  Technical  and 
engineering  support  within  areas  of  competence  is  provided  to 
ESD  Program  Offices  (POs)  and  other  ESD  elements  to 
perform  effective  acquisition  of  Cs/  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability /  maintainability  and  compatibility. 
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